Corrective lenses for high redshift super novae 
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Weak lensing of high-redshift Type la supernovae induces an external disper- 
sion in their observed standard candle brightnesses, comparable in magnitude 
to the intrinsic dispersion for redshifts z > 1. The same matter fluctuations 
£j , responsible for the magnification of distant supernovae also generate shear in the 

c3 . images of background galaxies. We investigate the possibility of using lensing 

shear maps constructed from galaxies surrounding the supernovae as a means of 
correcting the lensing-induced magnification dispersion. We find that a consid- 
erable fraction of the lensing dispersion derives from sub-arcminute scales, which 
are not probed by shear maps smoothed on arcminute scales. We thus find that 
weak lensing shear maps will be of only limited value in reducing the weak lensing 
magnification fluctuations of supernovae. 



1. Introduction 

High redshift Type la supernovae provide an excellent means of studying the expansion 
history of the universe (Riess et al. 1998; Perlmutter et al. 1999). It is estimated that 
the intrinsic dispersion in supernova luminosities can be calibrated to « 0.15 mag, and 
perhaps in the future to 0.1 mag (Kim et al. 2000), making them excellent standard candles. 
For supernovae at redshifts z < 1, this intrinsic dispersion sets the limiting accuracy with 
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which supernovae may be used to measure distances. For higher redshifts (z > 1), however, 
gravitational lensing by random fluctuations in the intervening matter distribution induces 
a dispersion in supernova brightness comparable to the intrinsic dispersion (Frieman 1997; 
Holz & Wald 1998; Holz 1998; Wang 1999), degrading their value as standard candles. These 
magnification fluctuations have zero mean, and so may be averaged away with sufficient 
numbers of supernovae (Wang 2000). However, the additional dispersion means that more 
supernovae are required than for low redshift samples to achieve a given signal to noise. 

It would be of great utility to determine the gravitational lensing magnification of each 
individual supernova. This would allow a correction of the observed brightnesses of the 
supernovae, and therefore improve their use as standard candles. Such a correction would 
be equivalent to obtaining a larger sample of supernovae, for free. In addition, measuring 
the gravitational lensing distribution at high redshift can be an important probe of the dark 
matter (Metcalf & Silk 1999; Seljak & Holz 1999; Metcalf 1999; Holz 2001). One means of 
achieving this would be an inspection of the foreground galaxies for each supernova. For 
example, SN1997ff at z = 1.7 has several foreground galaxies in its vicinity, leading to a 
magnification possibly as large as 0.4 magnitudes (Lewis & Ibata 2001; Riess et al. 2001; 
Moertsell, Gunnarsson, & Goobar 2001). If the magnification factor could be accurately 
estimated from the foreground galaxy images, then the supernova brightness could be cor- 
rected to its unlensed value. The correction factor depends strongly on uncertain properties 
of the galaxies' mass distributions (illustrated by the controversy over the extent of lensing of 
SN1997ff), and would miss possibly important contributions from dark halos. Furthermore, 
since such corrections would primarily shift highly magnified SNe to lower brightnesses, while 
leaving demagnified SNe unaffected, it would bias the resultant Hubble diagram. It is ap- 
parent that direct identification of individual lenses does not robustly determine the lensing 
magnification. It is also possible to correlate, in a statistical manner, the foreground galaxy 
number density close to the lines of sight to supernovae with the lensing effects on these 
supernovae (Metcalf 2001), but these statistical results do not help us "correct" any given 
individual supernova. 

An alternative method for correcting lensing magnification is to utilize weak lensing 
maps constructed from shear measurements of background galaxies. The same matter fluc- 
tuations responsible for the magnification of supernovae also lead to shearing of galaxy 
images. High redshift SNe are discovered by repeated exposures of wide fields, which when 
co-added provide extremely deep images of the galaxies surrounding the supernovae. Such 
deep, wide field images are well-suited for measurement of weak lensing shear. It is thus 
natural to hope that mass reconstruction from shear measurements of the surrounding fields 
might allow for the correction of weak lensing magnification, restoring the supernovae to their 
intrinsic brightnesses. A perfect measurement of the shear field at the redshift of a given 
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supernova would allow for a perfect reconstruction of the projected mass surface density 
(modulo the mass-sheet degeneracy, which should be unimportant for large enough fields). 
From this mass surface density it is possible to calculate the lensing magnification, and there- 
fore perfectly account for (and correct) the lensing effects on the observed brightness. Perfect 
shear maps are unavailable, however, and therefore our ability to infer the magnification is 
compromised. In this paper, we investigate how well weak lensing reconstruction can correct 
the brightnesses of distant supernovae. 

The basic scheme is as follows. A supernova occurs in a given field, and its peak apparent 
magnitude is observed and calibrated, using some variant of the Phillips (1993) relation. 
Then the (co-added) field containing the supernova is used to estimate the local shear at the 
supernova's location by averaging over a smoothing angle 9. The shear map is then converted 
to an effective convergence map using some reconstruction algorithm such as that of Kaiser 
& Squires (1993), and the derived convergence is used to correct the supernova's standard 
candle brightness. In the following section we estimate the variance in convergence for point 
sources given knowledge of the smoothed shear map, (k 2 ) 7 , which is a direct measure of 
the improvement such an approach can offer. We find that useful corrections require very 
large background source galaxy densities, and that this method is therefore of only marginal 
utility. 



2. Computing the shear-convergence correlation coefficient 

Let us denote by k the effective convergence, relative to the homogeneous filled-beam 
value, for a point source. In Figure 1 we plot the angular power spectrum of the convergence, 
Al(£) = £ 2 P K (£)/2ir (White & Hu 2000), for sources at z s = 2. To calculate this, we employ 
the fitting functions of Eisenstein & Hu (1999) for the linear matter power spectrum, and 
follow the prescription of Peacock & Dodds (1996) for the non-linear correction. We use a 
COBE normalized, scale invariant (n = 1) linear power spectrum in a flat ACDM cosmology 
with total matter density Q, m = 0.35, Hubble constant h = 0.65 (H = 100 /ikm/s/Mpc), 
and baryon density Vt b h 2 = 0.02. We also assume that the dark matter is microscopic (e.g., 
elementary particles), rather than macroscopic (e.g., black holes or MACHOs). The latter 
case leads to enhanced power on microarcsecond scales, which decorrelates point source 
magnification from galaxy shear. 

The convergence angular power spectrum peaks on arcminute scales (£ ~ 10 4 ), with 
significant power extending for multiple decades in i. All of this power contributes to the 
magnification of (almost) point sources like supernovae. When measuring shear, however, 
galaxy correlations must be averaged over large angular patches to suppress Poisson noise, 
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and this averaging washes out small scale power. For example, if we smooth over arcminute- 
sized patches, we see that a considerable fraction of the fluctuations affecting the brightness 
of point sources are not probed by the smoothed galaxy shear map. This hints that shear 
maps will be of only limited value. 

The convergence power spectrum gives the variance in effective convergence (White & 
Hu 2000) 

1 f°° 



4 V c 2 / Jo V a(R) J Jo k 2 

where R is the comoving radial distance (R(z) = dz' c/H(z')), a = 1/(1 + z), and 
^mass(^) a ) = k 3 -Pmass(^)/(2vr 2 ) is the matter power spectrum (per logarithmic interval phys- 
ical wavenumber). Here, and in what follows, we restrict ourselves to cosmologies with flat 
spatial sections (f2 tot = 1). For the power spectrum shown in Figure 1 we find (k 2 ) = 0.0036 
for sources at z s = 2. In the weak lensing limit, the magnification of a given source, /i, is 
related to the convergence by \x ~ 1 +2k. This variance in k thus corresponds to a la spread 
in standard candle flux of 12%. 

We define k 9 as the effective convergence averaged over a circular patch on the sky of 
radius 9, which is to be determined by shear measurements (e.g., using the algorithm of 
Kaiser & Squires (1993)). If (k 2 ) is the variance in k, then by measuring galaxy shear we 
can reduce the variance in the corrected convergence of the supernova image to (k 2 ) 7 = 
(1 — r 2 )(/t 2 ), where the correlation coefficient (in the absence of shot noise) is given by 
r 2 = (kkq) 2 / ((k 2 ) {kq)) . This expression strictly holds for a Gaussian convergence probability 
distribution function; although the weak lensing convergence pdf deviates somewhat from 
Gaussianity (see, e.g., Wang, Holz, & Munshi (2002)), we nonetheless expect the above 
expression for r 2 to be a reasonable ballpark estimate of the correction factor. It may be 
interesting to explore whether non-Gaussian pdf 's can lead to superior correction. Assuming 
we can write the shot noise term contribution to as Cp{9) = 7 2 ms /iV, with 7 rms the intrinsic 
galaxy ellipticity and iV = nir9 2 the number of source galaxies inside the smoothing area (n 
is the background galaxy number density), we have 

2 = (*Ke) 2 

~ {k 2 ){{k 2 ) + C p {9)Y 1 > 

For N oo and 9 — > 0, we have r 2 — > 1. This is as expected: with a perfect lensing shear 
map, we can fully correct for lensing magnification. 
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Fig. 1. — Weak lensing convergence angular power spectrum at z s = 2, for Q m = 1 — Q\ = 
0.35, h = 0.65, Q^h 2 = 0.02, and n — 1. We use the fitting functions of Eisenstein & 
Hu (1999) for the linear power spectrum (COBE normalized), and the fitting functions of 
Peacock & Dodds (1996) for the nonlinear power spectrum. Note that the convergence 
power spectrum peaks at angular scales of about an arcminute. Smoothing on these scales 
(or smaller) will therefore average away much of the convergence power. 
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The correlation functions can all be computed using the nonlinear matter power spec- 
trum (Jain & Seljak 1997; White & Hu 2000; Van Waerbeke et al. 2001), 

«)=T(¥^) 2 f SdR (SDTf^-*^^^ (4) 

and 

. , 9tt fVL m Hl\ 2 f Rs -R/R s )g(R) f°° dk 2 „ ,mwx,/, ™n ,r^ 

{KKe) = T{~^J I dR a{R? Jo ^ (*>«(*)) (5) 

where VT 2 (x) = 2J 1 (x)/x, J\(x) is the Bessel function of the first kind, and 

rRhor R — R 
g(R)=R dR s ^—w(R s ), (6) 

JR K s 

with w(R s ) describing the radial distribution of the sources and Rh OT the comoving distance 
to the horizon. If all of the source galaxies being utilized for the shear measurements are at 
the same redshift as that of the supernova (z s = zsn), we have w(R s ) = 5(R S — Rsn)- For 
a given background galaxy density, this case represents the optimal shear correction to the 
convergence. For the more realistic case of source galaxies distributed in redshift, we adopt a 
population distribution described by w(R s ) oc R° exp(— (R s /R+Y) (Kaiser 1992; Hu 1999), 
and set a — 1, (5 — 4, and R± = c/H corresponding to a mean redshift of z ~ 1. Although 
this galaxy distribution is simplistic (e.g., modulo an overall scaling factor it is independent 
of survey depth and filters), it is sufficient to indicate the decorrelation arising from the 
spread in galaxy redshifts. These two cases circumscribe the range of r 2 when tomographic 
information, such as photometric redshifts, are employed. 

For each value of the background galaxy density there is a tradeoff between the shot 
noise term, which decreases for large smoothing angles, and the cross-correlation term, which 
is weaker for large smoothing angles. The smoothing angle, 9, which maximizes the cross- 
correlation coefficient, r 2 , is shown as a function of n in Figure 2, for an intrinsic galaxy 
ellipticity of 7 rms = 0.4 (Kaiser 1992). The corresponding values of the cross- correlation 
coefficient are shown in Figure 3. The dashed curves represent the optimal correction, where 
all the source galaxies are at the same redshift as that of the supernova. The solid curves 
are for the galaxy redshift distribution described above. 



3. Conclusions 

From Figure 3 it is apparent that shear maps will be of limited value in reducing the 
lensing dispersion of supernova brightnesses, unless the number density of background galax- 
ies is great enough to permit sufficiently small smoothing angles. For example, if the number 
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Fig. 2. — The optimum smoothing angle, in arcminutes, as a function of the background 
galaxy density, chosen so as to maximize the correlation coefficient, r 2 . The supernova is at 
redshift z^n = 2, in a cosmology with Vt m = 1 — Qa = 0.35 and h = 0.65, and with an intrinsic 
galaxy shear of 7 rms = 0.4. The solid line employs the source galaxy distribution described 
in the text, while the dashed line represents the case where all of the source galaxies used to 
measure shear are at the same redshift as that of the supernova. 
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Fig. 3. — The correlation coefficient, r 2 , given by equation (3), as a function of the back- 
ground galaxy density. For each value of galaxy density the optimal smoothing angle is 
calculated (shown in Figure 2), and the subsequent shear correction to magnification is plot- 
ted. As in the previous Figure, the supernova is at redshift ^sn = 2, in a cosmology with 
Q m = 1 — Q\ = 0.35 and h = 0.65, and with an intrinsic galaxy ellipticity of 7 rms = 0.4. The 
solid line employs the source galaxy distribution described in the text, while the dashed line 
represents the case where all of the source galaxies used to measure shear are at the same 
redshift as that of the supernova. 
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density of background galaxies is as high as n = 10 6 /deg 2 , then at best we find a value for 
the cross-correlation coefficient of r 2 ~ 0.35. If the uncorrected convergence variance is 
(k 2 ) = 0.0036, then we find the variance for the corrected supernovae convergence to be 
(/t 2 ) 7 = (1 — r 2 )(/t 2 ) = 0.0023. This yields an rms magnification of 0.1, which is a 20% 
improvement over the uncorrected value of 0.12. This represents the optimal case discussed 
above, where all of the source galaxies are at the same redshift as that of the supernova. 
Including the expected spread in galaxy redshifts provides a more reasonable estimate of 
r 2 ~ 0.2, giving a reduction in the rms magnification of the supernova of around 10%. It is 
to be emphasized that this is an improvement (reduction) in the width of the observed su- 
pernova magnification distribution, and not a change in the mean (which remains at /i — 1). 
In addition, the intrinsic dispersion in supernova luminosities causes a further contribution 
to the observed rms standard candle magnification luminosity. 

Tomographic information can do little to ameliorate the situation. The simplest ap- 
proach would be to confine the shear analysis to source galaxies in a slab in redshift space 
centered on the supernova. By doing this one moves up and to the left of the solid curve 
in Figure 3, trading off increased shot noise for more effective lensing information. It is 
apparent from the Figure that there is a net improvement if > 10% of the galaxies are at 
similar redshifts to that of the supernova. More inspired schemes might attempt to employ 
the information contained in galaxies at all redshifts; regardless, for a given effective galaxy 
density the theoretical limit is still bounded by the dashed curve in Figure 3. Our model 
for the distribution of galaxy source images in redshift space is particularly simple — more 
realistic models (e.g. with dependence on survey depth) may also push one closer to the 
dashed curve. 

The estimate presented here is optimistic in that we assume that the smoothed con- 
vergence field may be directly measured. In reality, the shear field is measured, and then 
converted to a convergence map (Kaiser & Squires 1993). Even in this optimistic approxima- 
tion, at best meager returns are expected from the construction of shear maps of surrounding 
galaxies. An additional caveat is that we have assumed the noise is dominated by Poisson 
noise in the number of source galaxies. Additional systematic errors, such as imprecise 
measurement of the point spread function, only worsen the decorrelation. Note that our 
conclusions are sensitive to the shape of the convergence angular power spectrum. If in real- 
ity the power spectrum is unlike that of Figure 1, and instead has far less small-scale power, 
then galaxy shear may turn out to be a much more powerful tool for correcting weak lensing 
of supernovae. At present, it appears that there is significant small-scale power (Dalai & 
Kochanek 2002a,b), consistent with the values assumed here. Future wide field surveys like 
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the LSST 1 or SNAP 2 will directly measure the convergence angular power on some of the 
relevant scales, so it will be possible to check whether the assumptions made here are valid. 
Assuming, however, that the power spectrum does not significantly depart from that which 
we have used, the prospects for using galaxy shear to correct supernova brightnesses appear 
bleak. Given the danger of introducing unknown biases in the resulting distance-magnitude 
relation, it is unclear whether future supernova surveys should attempt the use of lensing 
shear maps to correct for magnification of supernova brightnesses. 
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